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■ Abstract. We review recent work on the Casimir interaction energy between 

cylindrical shells. We include proposals for future experiments involving cylinders, 



such us a null experiment using quasi-concentric cylinders, a cylinder in front a 
conducting plate, and a cylindrical version of the rack and pinion powered by Casimir 
lateral force. We also present an exact formula for the theoretical evaluation of the 
^ ' vacuum interaction energy between eccentric cylindrical shells, and describe improved 

^ , analytical and numerical evaluations for the particular case of concentric cylinders. 
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1. Introduction 

Up to now, most experiments aiming at a measurement of the Casimir force have been 
performed with parallel plates [1], or with a sphere in front of a plane [2]. The parallel 
plates configuration has a stronger signal, but the main experimental difficulty is to 
achieve parallelism between the plates. This problem is of course not present in the case 
of a sphere in front of a plane, but its drawback is that the force is several orders of 
magnitude smaller. On the theoretical side, the evaluation of the electromagnetic force 
for this configuration is still an open problem [3]. 

In this paper we will consider different configurations that involve cylindrical 
surfaces. As we will see, these configurations have both experimental and theoretical 
interest. In the next Section we will describe some promising experimental proposals. 
In Section 3 we will present an exact formula for the Casimir interaction energy between 
two eccentric cylindrical shells, and we will show that particular known formulas, such 
as the corresponding to a cylinder in front of a plane or the one of concentric cylinders, 
are included therein. Finally, for the particular case of concentric cylinders, in Section 
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Figure 1. Experimental setup for detecting Casimir forces using quasi-concentric 
cylinders. The inner cylinder is rigidly connected to a torsional balance and the signal 
to restore the concentric configuration is measured after a controlled displacement. 

4 we will present a new analytic result beyond the Proximity Force Approximation 
(PEA), and an improved numerical method to evaluate the interaction energy at small 
distances. 

2. Experimental proposals 

We discuss possible experimental arrangements for measuring the Casimir force between 
cylinders. As it was reported in the last years, cylindrical shells provide a new and 
promising arena to study Casimir interactions [H [5] . 

2.1. A null experiment 

Let us first consider the case of two eccentric cylindrical shells, in an almost coaxial 
configuration. The concentric is an unstable equilibrium position, so one possibility 
is to repeat a microscopic version of the experiment described in [6] to test universal 
gravitation in the cm range, with a small torsional balance mounted on the ends of the 
internal cylinder. In this case the unstable force could be evidenced by intentionally 
creating a controlled eccentricity and measuring the feedback force required to bring the 
internal cylinder to zero eccentricity, as depicted in Fig. 1 [1]. 

This configuration has some advantages over the parallel plates geometry. If there 
is no residual charge in the inner cylinder, the system remains neutral and screened by 
the external one from background noise sources, and from residual charges in the outer 
cylinder. When the inner cylinder has a residual charge, there will be a small potential 
difference between the cylinders, and the coaxial configuration will be electrostatically 
unstable. The electrostatic instability can be avoided by putting the cylinders in contact, 
something ineluctable during the preliminary stages of parallelization. Then the residual 
charge of the inner cylinder will flow to the hollow cylinder, apart from a residual charge 
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Figure 2. Cylindrical shell in front of a conducting plane. Using PFA one can show 
that the scaling of the Casimir force with distance is c>c d"^/^ [3], which is intermediate 
between the plane-spherical (cx d~^) and the parallel plate configuration (oc d"*). This 
configuration is also intermediate for the absolute value of the force signal for typical 
values of the geometrical parameters. 

due to the imperfections and finite length of the cyhnders. This residual charge will be 
smaller than for other geometries, as the same discharging procedure does not work in 
the other configurations (the efficiency of this procedure could be affected by the fact 
that repeated contact between the cylindrical shells may result in a degradation of the 
surfaces, as for instance an increased rugosity). 

The electrostatic instability could also be exploited to improve the parallelism 
between cylinders. One could apply a time-dependent potential between the cylinders 
and measure the force, as in the experiments to test the inverse-square Coulomb law. 
Parallelism and concentricity would be maximum for a minimum value of the force. 
Moreover, the expected gravitational force is obviously null, this being an advantage if 
one looks for intrinsically short-range extra-gravitational forces [7]. 

2.2. A cylinder in front of a conducting plane 

Another possibility, which in principle is much more appealing from the experimental 
point of view, is to consider the configuration of a cylinder in front of a plane. The 
cylinder-plane configuration of Fig. 2 is a compromise between the different drawbacks 
and advantages of the parallel-plates and the sphere-plate configurations. While this 
geometry offers a simpler way to control the parallelism with respect to the parallel 
plate case, at the same time gives rise to considerable force signal which requires the 
study of the force at larger distances. Not only can the study of the cylindrical-plane 
configuration provide insights into the thermal contribution to the Casimir force arising 
at any finite temperature, but also into the validity of the PFA. This is because, as the 
signal is stronger than in the sphere plane configuration, one can envisage sufficiently 
precise measurements at larger distances. This experiment is in progress [5], [8] . 

2.3. Cylindrical rack and pinion 

Another interesting configuration is a cylindrical version of the Casimir "rack and 
pinion" proposed in Ref.[9], in which the rack is a corrugated cylindrical shell that 
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Figure 3. Cylindrical version of the Casimir "rack and pinion" . Rotation of the outer 
cylinder induces a torque on the inner corrugated cylinder due to the lateral Casimir 
force. 



encloses the pinion, instead of a corrugated plane (see Fig. 3). In this case, the 
cylindrical geometry is of interest for the opposite reason than before: while reaching 
parallelism may be more complicated than from the plane rack, the torque will be 
enhanced by a geometric factor, and therefore could be more useful to generate the 
motion of the pinion. Indeed, the interaction energy per unit area between sinusoidally 
corrugated plane surfaces is given by [9] 



where d is the mean distance between the plates, h is the amplitude of the corrugations 
and X is the lateral displacement. We will not need the explicit form of the function J. 

The interaction energy for the plane and for the cylindrical rack and pinion can be 
easily computed using the PFA. In the first case it is given by 



where a and L are the radius and length of the cylinder, respectively, and d{d) = 
d + a(l — cos 9). 

For the cylindrical case we have, instead, E^rp = 2TTaLEpp. This simple result is 
valid when the radii of the outer (6) and the inner (a) cylinders satisfy a ~ 6 ^ d. 
Note that this configuration maximizes the superposition between the two surfaces, and 
therefore a uniform rotation of the external shell will produce a torque on the pinion 
much larger than the one produced by a plane rack moving with uniform velocity. 

The ratio of the forces can be easily estimated from the ratio of the interaction 
energies. Assuming that J(^y^) is a smooth function, and that the integral in Eq. ([2]) 
is dominated by 6* ~ 0, we obtain Fcrp/Fprp > ^fafd ^ 1. 




(1) 




(2) 
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Figure 4. Geometrical configuration for tlie eccentric cylinders. Two perfectly 
conducting cylinders of radii a < b, length L, and eccentricity e interacts via the 
Casimir force. 

3. The exact formula for eccentric cylinders 

The evaluation of the Casimir interaction energy between eccentric cylindrical shells 
(Fig. 4) has been initially performed using the PFA in Ref . |1] , where the force between 
a cylinder and a plate has also been computed and discussed for the first time using the 
same approximation. However, it is possible to go beyond the PFA, and find an exact 
formula for the interaction energy [TOl [TT] . This can be done using a mode by mode 
summation technique combined with the argument theorem. 

We start by expressing the Casimir energy as £" = {h/2) X]p(^p ~ ^p)^ where uip are 
the eigenfrequencies of the electromagnetic field satisfying perfect conductor boundary 
conditions on the cylindrical surfaces, and ujp are the corresponding ones to the reference 
vacuum (cylinders at infinite separation). In cylindrical coordinates, the eigenmodes are 
hn,k, = Rn{r,0)exp[-i{un,kj - hz)], where Up = Un,k, = v^/c^ + A^, and i?„, (A„) are 
the eigenfunctions (eigenvalues) of the 2D Helmholtz equation. Using the argument 
theorem the sum over eigenmodes can be written as an integral over the complex plane, 
with an exponential cutoff for regularization. In order to determine the part of the 
energy that depends on the separation between the two cylinders it is convenient to 
subtract the self-energies of the two isolated cylinders, -^12(0, b,e) = E — Ei{a) — Ei[b). 
Then the divergencies in E are cancelled out by those ones in Ei{a) and Ei{b), and the 
final result for the interaction energy is 

hcL f°° 

E^2{a,b,e) = ^ dy y log M{iy). (3) 

Here M = {E / E^o) / [{Fi{oo) / Ei{a)) (Fi(oo)/Fi(6))]. The function E is analytic and 
vanishes at all the eigenvalues A„ (-F00, at A„), and, similarly, Ei vanishes for all 
eigenvalues of the isolated cylinders. The function M is the ratio between a function 
corresponding to the actual geometrical configuration and the one with the conducting 
cylinders far away from each other. It is convenient to subtract a configuration of two 
cylinders with very large and very different radii, while keeping the same eccentricity of 
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the original configuration. Eq. ([3D is valid for two perfect conductors of any shape, as 
long as there is translational invariance along the z axis. 

The solution of the Helmholtz equation in the annulus region between eccentric 
cylinders has been considered in the framework of classical electrodynamics and fluid 
dynamics. The eigenfrequencies for Dirichlet boundary conditions (TM modes) and for 
Neumann boundary conditions (TE modes) are given by the zeros of the determinants 
of the non-diagonal matrices 

Q™ = [Jn{Xa)Nm{Xb) - Jm{Xb)Nn{Xa)] Jn^UXe), 
= [J:,{Xa)NUXb) - JiiXb)K{Xa)] Jn-m{Xe), 

where J„ and are Bessel functions of the first kind. The function M can be written 
as M = M^^M™, where M™ is built with 

= det [g™(a, 6, e)Q™{b, R, 0)] J] UXa), 

n 

F™{a) = det [Q™ia, R, 0)] J] ^(Aa), (4) 

n 

and i? is a very large radius. Similar expressions hold for M"""^. 

The Casimir energy can be decomposed as a sum of TE and TM contributions 



logM^^l^j+logM™ 



(5) 



with mte.tm(^) _ det[6np - A^™]. The non-diagonal matrices A™ and A™ are 

^np - 2^ J^(«;3) im-nmim-pm. 

Here J„ and Kn are modified Bessel functions of the first kind, a = b/a and 6 = e/a. 
The determinants are taken with respect to the integer indices n,p = — oo, . . . , oo, and 
the integer index m runs from — oo to oo. Eq.Q is the exact formula for the interaction 
Casimir energy between eccentric cylinders. 

This expression is rather complex to evaluate numerically, since each term in the 
infinite matrix is a series involving Bessel functions. However, we have been able to 
numerically evaluate the exact Casimir interaction energy between eccentric cylinders 
as a function of a for different values of the eccentricity |TT] . We have calculated this 
energy as a function of a for values that interpolate between the PEA (small a) values, 
and the asymptotic behavior for large a (see Eig. 5). The numerical convergence is 
better as a increases, while bigger matrices and more terms in each matrix element are 
needed as a gets closer to 1. 

An interesting property of the exact formula is that reproduces the exact Casimir 
interaction energy for the cylinder-plane configuration as a limiting case. Indeed, the 
eccentric cylinder configuration tends to the cylinder-plane configuration for large values 
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Figure 5. Exact Casimir interaction energy difference \AE\ between the eccentric 
and concentric configurations as a function of a ~ b/a for different values of S — e/a. 
Energies are measured in units of L/Aira?. 



of both the eccentricity e and the radius h of the outer cyhnder, keeping the radius a of 
the inner cyhnder and the distance d between the cyhnders fixed. Using the addition 
theorem and uniform expansions for Bessel functions it can be proved that, for x ^ h, 
Km{x + h) 



± rj 



{x + h) 



-Kn+p{2h). 



[x + h) 

Using these equations (with x = fie/ a and h = (3H/a) in our exact formula we find 



/iTM,c-p 

np 



^TE,c-p 



I'niP) 



Km 



Kn+p{2(3H/a), 
K,,+p{2f3H/a) 



which is the known result for the Casimir energy in the cylinder-plane configuration 



4. Concentric cylinders: new analytic and numerical results 

The exact formula for eccentric cylinders coincides, of course, with the known result for 
the Casimir energy for concentric cylinders (e = 0). Indeed, as /„_m(0) = 6nm, in this 
particular case the matrices A™'™ become diagonal and the exact formula reduces to 

T 

dp pluM'^'iP), (6) 



7-lCC 

^12 



where 



IniP)Kn{aP) 

h{aP)Km 



1 - 



I^{aP)Km 



(7) 
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The first factor corresponds to Diriclilet (TM) modes and the second one to Neumann 
(TE) modes. The concentric-cylinders configuration is interesting from a theoretical 
point of view, since it can be used to test analytic and numerical methods. It also has 
potential implications for the physics of nanotubes [T0| [T^. 

The short distance limit a — 1 <^ 1 has already been analyzed for this case [I3], and 
involves the summation over all values of n. As expected, the resulting value is equal 
to the one obtained via the proximity approximation, namely 

pec _ 1 /o^ 

%2,PFA- 3gQ^2 (c,_l)3- 

In the opposite limit (a ^ 1) it is easy to prove that to leading order only the n = 
term contributes to the interaction energy, and the energy decreases logarithmically with 
the ratio a = b/a, 

^ 1.26L 

^''^ Snb^lna 

It is worth noticing that, while for small values of a both TM and TE modes 
contribute with the same weight to the interaction energy, the TM modes dominate in 
the large a limit. 

4.I. Beyond proximity approximation: the next to next to leading order 

In this Section we will compute analytic corrections to the PEA given in Eq.Q. Due 
to the simplicity of this configuration, we will be able to obtain not only the next to 
leading order, but also the next to next to leading contribution. In order to do that, we 
need the uniform expansion of the Bessel functions. Eor example, we write 



,iy) ^ yrr^^-m ^^-^ ; ^ ^ ; ^ (n) 



where 



and similar expressions for the functions /.„. 

With these expansions at hand, we can evaluate the the matrix M both for the TE 
and TM modes. After a long calculation, it is possible to show that the Casimir energy, 
beyond the proximity approximation can be written as 

« -ssl^ ^ - 1' - - ^) (" - - ■■■} ■ (1^' 

In the expression above, the first term inside the parenthesis corresponds to the 
proximity approximation contribution in Eq.Q, while the second and third terms are 
the first and second order corrections, respectively. It is important to stress here that 
both TM and TE modes contribute with the same weight to the energy up to the next 
to leading order, but it is not the case in the third term. There is a factor I/tt^ coming 
from the TM mode, and a factor I/tt^ + 1/10 corresponding to the TE one. 
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Figure 6. Ratio between tlie exact Casimir energy for concentric cylinders and 
the Casimir energy estimated using the PFA up to the next to leading order Epf^, as 
a function of the parameter a. This is done for two different methods: the numerical 
(of slow convergence) and the numerical improved (subtraction method). Solid line 
indicates the analytic result from Eq. fT^ . 



4-2. Improving the convergence of the numerical evaluation 

Numerical calculations of the Casimir energy for a very close to one are really difficult 
since big number of terms have to be considered in the sums, and therefore convergence 
problems arise. 

In order to perform a numerical evaluation of the Casimir energy for the concentric 
cylinder case, we will describe a subtraction method, in which we have used the proximity 
approximation value of the energy to improve numerics. 

As a guiding example, let us consider the following sum: 

M ^ 
n=l 

The convergence of this sum as M — oo is extremely slow, as the following numbers 
suggest: ziQS = 5.5728, ziqs = 7.4222 and z^o = 10.5844. About 10^° terms are needed 
to get an accuracy of 1%. 

Then, we add and subtract the function J^^ p^, that reproduces the behaviour of 
the series when the exponent in the denominator is close to one. In this way, 

f'' dx dx f'^ dx 

ZM = ZM-j^ ^1=Dm + J^ ^-^oo + lO. (14) 

The convergence of the new function Dm is notably faster as Dio = 0.6234 and 
-Diooo = 0.5847, i.e. 1% accuracy is obtained with less than 10 terms. 

In the case we are concerned here, we can add and subtract the interaction energy 
for concentric cylinders computed using the leading uniform asymptotic expansion for 
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the Bessel functions, up to first order in a — 1: 

K^jnay) Injny) ^ ^_2nia~i)y/Iw _ (I5) 
Kn{ny) In{nay) 

We denote by E the interaction energy obtained by inserting these expansions into Eq. 
([6]). Now we write 

Ell = {Ell -E) + E. (16) 

The last term in the above expression is easily written into an analytic expression, which 
contains the leading order of the Casimir energy. Meanwhile, the difference contained in 
the brackets in Eq. (fT6!) . has a faster convergence than the original sum and therefore, 
can be easily calculated numerically. 

In this context, in Figl6]we present both Casimir energy of the concentric cylinders 
for the direct numerical calculation (of slow convergence) and the alternative method 
mentioned above. In this figure we plot the ratio Ei2/Epf^ where 

As can be seen, with this subtraction method it is possible to compute the exact energy 
for values of a much closer to 1, while the accuracy of the direct calculation is worse for 
a < 1.02 . 

A similar method could in principle be applied to the eccentric cylinders or the 
cylinder-plane configurations, although in these cases the main difficulty is the analytic 
evaluation of the approximate energy to be added and subtracted. 



5. Final remarks 



In this paper we have described experimental proposals and theoretical aspects of the 
Casimir interaction energy between cylindrical shells. We have reviewed previous works 
on the subject, in which we obtained an exact formula for the interaction energy between 
eccentric cylinders (TUl [H], and where we discussed the advantages of considering 
experiments involving cylinders [ll[5]. 

We have also presented some new results. In particular, we have shown that 
a cylindrical version of the non contact rack and pinion powered by the lateral 
Casimir force, proposed in Ref. [9], would have a larger torque because of a geometric 
enhancement. 

From a theoretical point of view, we have found an analytic formula for the 
interaction energy between concentric cylinders beyond the PEA, including first and 
second order corrections. We have also presented a subtraction method useful to improve 
the convergence of the numerical calculations as the concentric surfaces get closer to each 
other. We hope to generalize these results to other geometries in future works. 
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